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SUMMARY 

Hie NASTHAN THAPRG and THAPAX finite elanents are very restrictive as to shape 
and grid point numbering. Hie elements must bo trapezoidal with two sides parallel 
to the radial axis. In addition, the ordering of the grid points on the element 
connection card must follow strict rules. Hie paper describes the generalization of 
ttiese elements so that these restrictions no longer apply. 


INTRODUCTION 

Since NASTRAN's inception in the early 1970's, the axisymmetric trapezoidal 
ring element TRAPRG has been an accurate, efficient element used in solid arisymmet- 
ric problems. More recently, the THAPAX elenent was introduced to handle non- 
axisymmetric loading for such structures. While these two elements usually perform 
very well, the restrictions imposed upon the user in the specification of the 
elements can be considered, at best, difficult, at worst, unreasonable, in light of 
today's automatic data generators. Hie restrictions require that the elements be 
trapezoidal, rather then genercily quadrilateral, with the top and bottom edges 
parallel to the radial axis. Also, the specification of the grid points on the 
connection card i ust be given counterclockwise starting with the grid point with 
smaller radial coordinate of the two grid points with the smaller axial coordinate. 
Hiis paper presents the reasons for these restrictions, how the restrictions have 
been removed, the imposition of a new, out less stringent restriction, and an 
exanple problem. 


THEORY 

Hie stiffness matrix for a finite elenent is usually represented as 

[K] = J v [Bl T [D][B]dv (1) 


where LB] is the matrix of strain-displacement relations and [D] is the materials 
matrix describing the constitutive relations. The NASTHAN Theoretical Manual (ref. 
1) defines LB] for THAPRG and THAPAX elements. The manual shows that, in order to 
evaluate the integral of Equation (1), the integrals 


*A r P z q drdz, p = -1 , 0, 1 , 2, 3 and q = 0, 1 , 2 (2) 
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must be evaluated over the cross-sectional area of the element. The NASTRAN Pro- 
grammer's Manual (ref. 2) describes how the integrals (2) are computed analytically. 
It is because of the nature of the expressions involved in the exact representation 
of these integrals that the original restrictions on the elements were imposed. The 
restrictions have been removed by replacing the exact integration with Uauss quadra- 
ture, as follows. 


Ihe three-point quadrature is given as 
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f k r*VW . £ ^ » ) H,H U(5 1 , n >1 
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N. = 1/4(1 ♦ £5.)(1 + nn i ) , the linear isoparametric shape function 
over the square £ = -1 to +1 and r» = -1 to +1 


(3) 




r^ = r.z coordinates of the grid points at the four corners of 
the elaaent 
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With this formulation, the restrictions on the trapezoidal shape and grid point num- 
bering can be ranoved. 

There is, however, one situation which can cause numerical problems. Note 
that, in integrals (2), if p = -1 and r -*■ 0, then the integral -*■ “ . Two cases can 
be examined. The first is the core element, a TRAPRG element with exactly two grid 
points with zero radius. (NASTRAN does not allow the TRAPAX elaaent to be a core 
element since it does not allow a RINUAX ring to have a zero radius.) At these two 
| points, the radial displacaoent must be zero. This condition can be assured as 
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follows. If u and w are the radial and axial displacements, respectively, at some 
point within an element, then the assumed displaceaent functions are (ref. 1) 


U = B 1 

♦ &2 r * 

♦ Vz 

(A) 

W =B 5 

+ Qy-r + B„z 
o 7 

♦ B a rz 

(5) 


where the are unknown coefficients. Ihen, 



where u. , w. are the displacements of the ith grid point. If r. = 0, then u. must 
be 0, which x means that the (2 -1)th col min of H can be zeroed x In addition, from 
Equation (4), 0 = u s B ♦ B z for all z. Therefore, 0 = B. = B , which can be 
assured by zeroing the firsthand third columns of H . However, the only terms of 
hL which contain a form of (2) with p = -1 are in uie first and third columns. 
Therefore, for a TRAPRG core element with r, = 0, no nuaerical problaos exist. 


Now consider a second case - one where some r^ is small compared with some 
other r. in the elaaent. Comparisons between the analytical integration and the 
Gauss quadrature show that if 


IsiJi > io (7) 

r inin 

then the Gauss quadrature results begin to significantly degrade (r = max r. , 
i = 1, ..., 4). Therefore, inequality (7) is not allowed. This restriction should 
impose a hardship mainly when a TRAPAX core element is desired. This hardship how- 
ever should require only an extra eleoent or two to transition to larger radii. 


SAMPLE PROBLEM 

The sample problaa for this work was a normal modes analysis of a thick-walled 
cylinder of inner radius 5.0 inches, outer radius 6.25 inches, and length 3.5 
inches. The finite elaaent .model for the unaodified version of NASTRAN is shown in 
Figure la and, for the version which removes the restriction, in Figure 1b. Both 
models used the TRAPAX element. The results are shown in Table 1. The comparison 
of the results shows a degradation at the highest mode computed . However, this can 
be expected because of the severity of the non-uniform mesh. 
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TABLE 1. COMPARATIVE RESULTS 


NATURAL FREQUENCIES (CPS) 


MODE 

RESTRICTIVE 

ELEMENTS 

GENERAL 

ELEMENTS 

1 

1090. 002 

1090.000 

2 

1660.250 

1660 . 864 

2 

1760. 024 

1760. 101 

4 

1968. 622 

1989. 121 

5 

2071.862 

2071. 920 

6 

5407.211 

5209.751 
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